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ABSTRACT 

High-energy cosmic rays passing through the Earth's atmosphere produce extensive showers whose charges emit radio 
frequency pulses. Despite the low density of the Earth's atmosphere, this emission should be affected by the air refractive 
index because the bulk of the shower particles move roughly at the speed of radio waves, so that the retarded altitude 
of emission, the relativistic boost and the emission pattern are modified. We consider in this paper the contribution of 
the boosted Coulomb and the Cerenkov fields and calculate analytically the spectrum using a very simplified model in 
order to highlight the main properties. We find that typically the lower half of the shower charge energy distribution 
produces a boosted Coulomb field, of amplitude comparable to the levels measured and to those calculated previously 
for synchrotron emission. Higher energy particles produce instead a Cerenkov-like field, whose amplitude may be smaller 
because both the negative charge excess and the separation between charges of opposite signs are small at these energies. 

Key words, radiation mechanisms: non-thermal - waves - elementary particles - cosmic rays 



1. Introduction 

Nearly a century after the discovery of cosmic rays, their 
nature and origin still constitute one of the great problem 
of contemporary astrophysics. Most of them are revealed by 
their electromagnetic radiation, but those reaching the solar 
system play a unique role since they can be observed nearly 
directly. The extensive air showers made of secondary (or 
higher order) particles produced by high-energy cosmic rays 
in the Earth's atmosphere are currently studied with giant 
ground-based particle detectors, whereas additional infor- 
mation is provided by observing fluorescence and Cerenkov 
radiation in the optical range. 

A complementary technique is being developed, based 
on the detection of radio emission from the shower charged 
particles, making advantage of the techniques of increas- 
ing sophistication developed for radioastronomy. Indeed, it 
has been known for several decades that the charges in ex- 
tensive cosmic ray showers produce radio frequency pulses, 
as first suggested by Askar'yan (|1962p who estimated the 
Cerenkov emission. Subsequent studies (see Kahn & Lerche 
I1966[) have suggested radio emission to be produced instead 
by acceleration by the Lorentz force in the Earth's magnetic 
field (see the extensive review by Allan I197ip . and have 
stimulated extensive calculations and modelling of the syn- 
chrotron radiation (see Huege & al. 120071 and refs. therein). 
Up to now, the difficulties in both theory and observation 
have precluded a full understanding of the origin of the 
radio emission. 
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Most previous studies of this radio emission have ap- 
proximated the air refractive index by unity. However, the 
median energy of electrons in the shower is about 30 MeV, 
corresponding to 1 — v/c ~ 1.5 x 10 -4 ; the median speed 
is thus greater than the phase speed c/n of radio waves in 
air at sea level (where n ~ 1.0003), and is roughly equal to 
the phase speed at an altitude of half the atmospheric scale 
height. The air refractive index may thus affect significantly 
the radio electric field produced by the shower, whatever 
the emission mechanism, by changing the retarded altitude, 
the relativistic boost and the radiation pattern. 

In particular, this should increase the relativistic boost 
of the Coulomb field produced by the shower charges, at 
speeds slightly smaller than the wave phase speed, whereas 
faster charges, which can catch up with the waves they emit, 
produce a Cerenkov like field. We make below an analytical 
calculation of these fields with a highly simplified model of 
the charge distribution, in order to highlight the physical 
processes and evaluate the contributions to the total electric 
field. 

Before taking the refractive index into account, let us 
estimate the boosted Coulomb field produced in vacuum by 
a shower particle. Consider a charge q in uniform relativistic 
motion. For an observer at rest, the Coulomb field is radial 
relative to the instantaneous present position of the charge, 
but the FitzGerald-Lorentz contraction compresses the field 
lines, so that in the directions normal to the velocity, the 
field in vacuum is greater than the isotropic Coulomb field 
by the Lorentz factor 7. This produces a pulse of elec- 
tric field of maximum amplitude Ej^ ax = ^q/ Aire^d 2 when 
the charge's path passes at closest distance d to the an- 
tenna, and of temporal width r = d/^c), directed from 
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the charge's present position to the observation point. The 
number of charges in the shower at maximum develop- 
ment is about one per GeV of primary energy with about 
20% relative excess of electrons over positrons. With a pri- 
mary cosmic ray of 2 x 10 17 eV, i.e. 2 x 10 8 charges, and 
a median secondary energy of 30 MeV (7 = 60), this 
yields a pulse of maximum amplitude a few 10~ 4 V/m 
and of duration roughly 10 ns at 100 m perpendicular dis- 
tance. The corresponding low-frequency Fourier transform 
ZtEmcix — 9/(27reorfc) is a few juV/m/MHz at frequencies 
smaller than l/27rr ~ 50 MHz. The exact strength and po- 
larisation are determined by the magnetic charge separation 
which produces oppositely directed changes of position and 
velocity for charges of opposite signs. This field should be 
added to the contribution of synchrotron emission. 

Since the above amplitude is of the same order of mag- 
nitude as the values observed by different instruments (see 
Allan [19711 Ardouin [20061 Horneffer [20061) , this contribu- 
tion may not be negligible. We therefore consider this field 
in more detail, taking the refractive index into account. 



2. Electric field spectrum of a subluminal charge 

In this Section, we consider a point charge. A charge dis- 
tribution and the corresponding coherence effects will be 
considered in Section 4. 



2.1. Retarded potentials 

Consider a point charge q moving in a medium of constant 
refractive index n (and relative magnetic permeability /1 ~ 
1), at velocity v. The field potentials at space-time point 
(r, t) are the standard Lienard-Wiechert potentials with the 
formal replacements c — > c/n and eo — > eon 2 



*(r,t) 



1 



1 - n/3 ■ n I R 



A(r,t) = n 2 /3$(r, t)/c 



(1) 
(2) 



where P — v/c, n is the unit vector in the direction of 
the point of observation from the moving charge, R is the 
distance to the moving charge, and the subscript "ret" on 
the square bracket indicates that the quantities inside are 
to be evaluated at the (observer's) retarded time t ret =t — 
nR(t ret )/c (see for example Fevnman 119641 Jackson [1999, 
ThideQMZD- 

Equations (HJ-© hold whatever the charge's motion 
and the refractive index, provided n is a constant. Namely, 
this formulation neglects both the variation of n with po- 
sition and the dispersion, in particular the variation of n 
with frequency (Ginzburg 1989, Clemmow & Dougherty 



1994). In this paper we will apply these equations for a 
uniform velocity, so that the electric field, given by (|15[) . is 
the acceleration-independent term of the Lienard-Wiechert 
electric field (see for example Fevnman ll9641 Jackson[1999), 
i.e. mainly a boosted Coulomb field for a subluminal charge 
or a Cerenkov-like field for a supraluminal charge. This field 
should be added to the synchrotron field calculated by the 
current models of shower emission. In the present Section, 
we consider the first case: v smaller than the velocity of 
light in the medium c/n, i.e. n(3 < 1. 

In this case, only one instant in the charge's past history 
has a light cone that reaches a given location in space-time 
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Fig. 1. Present (P) and retarded (P') positions of the mov- 
ing charge with respect to the antenna (A) for n{3 < 1 (not 
to scale). 



(r, t), so that there is only one retarded time t re t for a given 
(r,i). Let us take the origin of (observer's) co-ordinates so 
that the charge moving along the z axis passes at z = at 
t = 0, and let us put an electric antenna at co-ordinate z 
along the charge's path, and distance d in the perpendicular 
plane along the x axis (FigfT]) . This choice of origin allows 
to keep the time of arrival of the charge in the results, which 
will be used in Sections 4 and 6 to evaluate coherence effects 
for a source of finite size. 

The distance between the charge's present and past po- 
sitions is PP'= n(3R, so that with the notations of Fig[Tl 
P'Q= n (P ■ n) R, whence AQ= (1 - n/3 • n) R. From the 
triangles CPA and QPA we have d 2 + (z + vt) 2 =AQ 2 + 
(n(3Rsm9) 2 , so that, since sinf5 = d/R, the square of the 
denominator in the bracket of (fl]) is 



[(1 - np • n) R] 2 ret = d 2 (1 - n 2 f) + (z + vt) 2 



(3) 



Equation |3]) is a purely geometrical relation which holds 
whatever the sign of (l — n 2 /3 2 ) , provided at least one point 
P' (source's retarded position) does exist. This is the case 
when the right-hand side term of ([3]) is positive or zero, a 
condition which always holds for the subluminal case (n/3 < 
1), and which defines the Cerenkov cone for n(3 > 1. We 
shall consider the latter case in Section 5. 
Substituting © into (JTJ) yields 



$(r,t) 



where 



47re n 2 d(l -n 2 f3 2 ) 



d\ l-n 2 /? 2 I 1 / 2 /v 



t + z/v 



-1/2 



(4) 



(5) 



Note that we have used an absolute value in (0, although 
1 — n 2 (3 2 > for a subluminal charge, in order to be able 
to use the same definition of r for a supraluminal charge 
(Section 5). 

Equations (g]) and © yield a pulse in the potentials, 
centred at t = —z/v, the time when the charge's trajectory 
passes at closest approach to the antenna, and of half- width 
about r. In the limit n — > 1, (HJ) reduces to the well-known 
potential of a point charge in uniform motion in vacuum 
(see for example Feynman[1964, Jackson [1999]). 
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It is important to recall that in ([?]), the relevant charge 
and velocity are retarded quantities. Namely, the potentials 
only depend on the charge, position and speed of the parti- 
cle at the retarded time. This means that the pulse's peak 
is produced by the past state of the charge, at the retarded 
time t ret , when the charge was at distance z re t — z = dj tan 9 
from the antenna along the direction of motion. Since at 
the time of the peak of potential at the antenna (A), the 
charge's position (P) is at C so that tanfJ = dj (nj3R), we 
have cos 9 — n(3 at this time, whence 



Zret 



i(3d (1 - n 2 f3 2 ) 



-1/2 



(6) 



The state of the charge at times different from the re- 
tarded time that produce the pulse is irrelevant. In partic- 
ular, the instant of the peak is that at which the charge 
moving at v along the z axis passes at closest approach 
to the antenna, but the real charge need not necessarily be 
there at this time; it may have disappeared or have changed 
its trajectory; the latter effects do not change the poten- 
tials, but produce an additional contribution to the electric 
field (dnj) due to the term dA/dt (see Feynman & al. H964p . 

For cosmic showers, therefore, the relevant charge and 
speed producing the field are those corresponding to the 
development of the shower at the retarded altitude, given by 
([6]); hence they do not necessarily correspond to maximum 
shower development. For an antenna on the ground and 
a shower inclined to the vertical by the angle a, we see 
from FigQ]and Eq.([6]) that the retarded altitude is h ret = 
GP'cosa, i.e. 



d n(3 (l — n 2 /? 2 ) cos a ± sin a 



(J) 



where the signs + and — correspond respectively to the 
charge impacting the ground after or before the time of 
closest approach. For 7 > 1, u - 1 < 1, CP'= z ret — z is 
much greater than GC, so that ([7]) yields approximately 



h ret ~ d (1 - n 2 (3 2 ) 



-1/2 



cos a 



(8) 



where we have also omitted the factor n(3 ~ 1. We shall 
return to the retarded altitude later. 

From ([?]), we deduce immediately the x component of 
the electric field E x — —d<&/dd 



4ire n 2 d 2 (1 - n 2 /3 2 



a/2 



1 + 



t + z/v 



-3/2 



(9) 



2.2. Effect of the refractive index 

For relativistic velocities, the refractive index plays an im- 
portant role, even though it is close to 1. Indeed, we see 

from HI)-© that a subluminal charge q of Lorentz factor 

r /2 

7 = (l — /3 2 ) in a medium of refractive index n ^ 1 
produces the same electric potential as a charge q/n 2 of 
Lorentz factor 



7 



2 «2\ 



1 - n z (3 



-1/2 



(10) 



in vacuum. In particular, the pulse maximum amplitude 
and retarded altitude are both proportional to 7', whereas 
the time duration is inversely proportional to 7'. 

Figure [2] shows the factor 7' as a function of 7 for dif- 
ferent values of the air refractive index. One sees that, 




n = 1.0002 n = £.0001 



10 20 30 40 50 60 70 80 
Y 

Fig. 2. A subluminal charge of Lorentz factor 7 in a 
medium of refractive index n produces a similar electric po- 
tential as a charge in vacuum of Lorentz factor 7', plotted 
here versus 7 for different values of the air refractive index. 
When 7 approaches the value for which n(3 =1,7' increases 
rapidly and both the variation of n and the strong increase 
in retarded altitude come into play. For still greater values 
of 7, 7' is imaginary and Cerenkov radiation is produced. 



as we noted in the introduction, the subluminal upper 
limit nj3 = 1 takes place at an energy of the same or- 
der of magnitude as the median energy of shower charges 
(7 ~ 60). Note that for 7 3> 1, n — 1 -C 1, we have 
1 — n 2 (3 2 ~ I/7 2 — 2 (n — 1), so that 7' may be approxi- 
mated by 



7' ~ 7 [1 - 2 (n - 1) 7 2 



-1/2 



(11) 



In practice, since the refractive index decreases as al- 
titude increases, the retarded distance R should be eval- 
uated by an integration involving the wave phase speed 
at different altitudes from emission to reception. Since the 
greater the refractive index, the greater the propagation 
time, this integral is mainly determined by the greater val- 
ues of the index. To simplify the calculations, we thus ap- 
proximate the index by its value at low to mid altitudes, 
n ~ 1 + 2 x 10~ 4 , but will keep in mind that n — 1 is 
roughly proportional to the air density, decreasing roughly 
exponentially with altitude, and increasing with humidity 
(see for example Birch I1994p , which may have important 
consequences. 

With this value of n, we see from (fTTj) that 7' ~ 7 only 



for 7 significantly smaller than [2 (n — 1)] 



-1/2 



50, i.e. 



for electrons of energy significantly smaller than about 25 
MeV. Since the median electron energy in a typical shower 
is about 30 MeV (7 ~ 60), the radio electric field is expected 
to be rather different from its value in vacuum. For example, 
electrons of 7 ~ 35 (respectively 40) in air with the above 
value of n produce a similar electric field as electrons having 



7 



49 (respectively 67) in vacuum. 



Close to n/3 — 1, i.e. 7 ~ [2 (n — 1)] , 7' varies very 
rapidly. In this region, two physical effects come into play: 

— the variation of the refractive index, since 7' is very 
sensitive to n in this region, 

— the strong increase of the retarded altitude with 7', 
which tends to put it above the atmosphere, so that 
the retarded charge becomes negligible. 

For example, at perpendicular distance d = 200 m and 
vertical inclination angle a = ir/4, ([5]) yields h ret = 20 km 
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for 7' ~ 140, i.e. 7 ~ 47. Hence in that case, for the retarded 
altitude not to be above the region of shower development, 
the energy of the radiating electrons should be smaller than 
about 24 MeV. This means that, since the median energy 
is about 30 MeV), typically a little less than the lower half 
of the electron energy distribution in the shower contribute 
to the boosted Coulomb electric field, for these values of 
perpendicular distance and angle to the vertical. 

For greater energies, the retarded altitude is above the 
atmosphere, except for showers of large inclinations or pass- 
ing very close to the antenna. And at energies so that 
n/3 > 1, 7' becomes imaginary; in that case, there is no 
longer a single retarded time and position, and a Cerenkov 
field is produced, which we shall evaluate in Section [5J 

2.3. Electric field spectrum 

Defining the Fourier transforms of the potentials as 

f + OO 



/TOO 
dte luJt <f>(r,t) 
-00 

we have from (J4j and ([2]) 



*(r,w) 



2ireon 2 v 
A(r,w) = n 2 (3$(r,u)/c 



K {ud/^'v) 



(12) 

(13) 
(14) 



Here Kq is a modified Bessel function of order 1 
(Abramowitz & Stegun fl972|l . r is given by (0, z is the 
antenna's co-ordinate along the charge's path (whose origin 
is the charge's co-ordinate at t = 0), and d is the antenna's 
perpendicular distance to the charge's path. 
The electric field is given by 



E(r,<) = - xj<!>-dA/dt 



(15) 



so that its Fourier transform has from (fl"3"|) - ([T4"|) the com- 
ponents 



E x (r,u) 



E z {y,u) 



que 



z jv 



V ' 7 



-ique 



iujz/v 



27reon 2 i' 2 7' 2 



K Q {ud/iv) 



(16) 
(17) 



Therefore E z /E x <lfor7>l,n-l<l, so that the 
electric field is radial (perpendicular to v, directed along 
the charge's present position to the antenna) of amplitude 



E(r,u) 



que 



iuz/v 



7 Ki {ud/iv) 



(18) 



with 7' given by (|10[) . at perpendicular distance d from the 
charge's path. 

Using the expansions of the Bessel function K\ 
(Abramowitz & Stegun I1972p . (fT5|) yields at respectively 
low and high frequencies (or distances) 



ud , qe l " z l v 

— « 1 E(t,u) ~ f — 

7'u Zireon^vd 

ud „, , qe iwz / v /2nui 

— » 1 E(r,u) ~ f — —— 



1/2 



odf-y'v 



(19) 



(20) 



for n— l<Cl,/3 — 1. The low- frequency (or small distance) 
spectrum (the time integral of the pulse) is independent of 
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Fig. 3. Spectral field strength | E(u) \ from (JT5J), produced 
by a point charge of 10 8 electrons at perpendicular dis- 
tances d = 100 and 200 m, for n = 1.0002 and 7 = 40. 

7', as expected since the larger 7', the greater the amplitude 
of the pulse, but the smaller (in the same proportion) the 
duration. At large frequencies and/or distances, the field 
decreases nearly exponentially with the product ud/^'v, 
thus with a frequency scale proportional to 1 j d and a dis- 
tance scale proportional to 1/u. At intermediate frequen- 
cies, the factor \Jujd in (|2H)) makes the decrease with u 
(respectively d) slower (respectively faster) than given by 
the exponential e~ ud ^ v alone. 

Figure [3] shows the modulus of the electric field spec- 
trum (fTSf produced by a charge q = Ne with N — 10 8 
(e being the electron charge) for n = 1.0002 and 7 — 40 
(7' ~ 67), at perpendicular distances d = 100 and 200 m, 
for which the retarded altitude ([8]) is respectively 5 and 10 
km at a vertical angle a — tt/4. 

Note that the plotted values correspond to Fourier 
transforms (defined as in (fT2)) ), so that they should be mul- 
tiplied by a factor of two when comparing them to measured 
spectra, which are generally defined for positive frequencies 
only. 

3. Charge at the retarded altitude 

The retarded charge producing the pulse is that correspond- 
ing to the shower development at the retarded altitude, 
given by {8J as a function of the perpendicular distance 
d to the path and of the inclination a to the vertical. In 
the region of maximum shower development, the number 
of charged particles (mainly electrons and positrons) is 



N M 0.31 (W p /W c ) [In {W p /W c )\ 1/2 



(21) 



for a primary of energy W p , where W c = 0.86 x 10 8 eV 
(see for example Abu-Zayyad, & al. 120011 and references 
therein). 

With a simplified atmosphere model of density p = 
poe~ h ' H at altitude h, the total mass from the top of the 
atmosphere down to altitude h is 



X{h) = (p H/cosa)e- h/H 



(22) 



where H ~ 8500 m and p ~ 1.22 kg/m 3 . The altitude 
Km where the shower charge is maximum is determined 
by X(Iim) = Xm, where Xm — 6500 kg/m 2 for a primary 
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Fig. 4. Number of secondary charges in the shower (nor- 
malised to the value at maximum shower development) at 
the retarded altitude 7'dcosa for 7 = 35 (solid lines) and 
7 = 40 (dotted lines) with n — 1.0002 (i.e. respectively 
7' = 49 and 67) as a function of distance for different an- 
gles a to the vertical, for a primary of 10 17 — 10 18 eV. 



cosmic ray of energy 10 17 — 10 18 eV, and increases at higher 
energies (Abu-Zayyad, & al. 120011 and references therein). 
The number of charged particles at the retarded altitude is 
determined by the value of X. It can be obtained from the 
Greisen function as 



N(X)/N M = exp 



L C V 2 J L c 



(23) 



where the so-called shower age is s = 3/ (1 + 2Xm/X), Nm 
is given by (HU, and L c ^ 587 kg/m 2 ( Abu-Zayyad, & al. 
I20UT1) . 

Figure[4]shows N/Nm as a function of the perpendicular 
distance d for different values of the inclination a of the 
shower to the vertical, for 7 = 35 and 40 with n — 1 = 
2 x 10~ 4 (respectively 7' = 49 and 67). 

At normal incidence, the number of charges at the re- 
tarded altitude is of the same order of magnitude as the 
value at maximum development, up to distances of a few 
hundred metres, for which the retarded altitude goes above 
the region of significant shower development. As the angle a 
increases, the distance d for which the retarded altitude is in 
the region of maximum development increases. The larger 
7 and a, the larger the value of d for which the retarded 
altitude corresponds to maximum shower development. As 
7 increases significantly above 45 (for n — 1 ~ 2 x 10 -4 ), 
the fast increase of 7' puts the retarded altitude above the 
atmosphere. 

In the following Section, we will thus consider charges 
of Lorentz factor 7 = 40. 

4. Boosted Coulomb spectrum of a charge 
distribution 

In the previous Sections, we have calculated the electric 
field of a point charge. However, the charge distribution in 
the shower is believed to be pancake shaped, of axis the 
shower direction, and of lateral size increasing with alti- 
tude. The half-size containing the bulk of the particles is 
about 20-30 m at ground level, somewhat smaller than the 
canonical Moliere radius of about 80 m (Antoni & al.[2001), 



d 



Fig. 5. Geometry for calculating the field at perpendicular 
distance d from the centre of a charge distribution with 
cylindrical symmetry around the charge path. 

and very thin in the longitudinal direction, with a disper- 
sion in particle arrival times of about 1.6 ns at the centre 
(Linsley 1986). This small longitudinal dispersion will de- 
crease the field level because of the loss of coherence due to 
the factor e luz '" in the spectrum (fT5|) . 

With an axial symmetry, it is convenient to consider 
first a charge distributed round a circle of radius a in the 
plane perpendicular to the path, infinitely thin in the z 
direction, which allows an analytical calculation. The field 
of an arbitrary charge distribution can then be obtained by 
performing an integration over the lateral and longitudinal 
dimensions. 



4.1. Ring charge distribution 

Consider a charge q distributed round a circle of radius a 
in the plane perpendicular to z which passes at z = at 
t = 0. The potential produced at perpendicular distance d 
from the centre of the circle (and co-ordinate z along the 
charge's path) is, from (fT3"|) with the geometry defined in 
FigE] 

a iujz I v 

-F(a) (24) 



$a(r,w) 



qe" 



27reon 2 w 
2n 



a 



d(f)Ko(ujr/"f'v) 

1 /2 

2 - 2ad cos 6) 



(25) 
(26) 

The integral F (a) can be calculated analytically by ex- 
panding Kq in series of Bessel functions using Neumann's 
addition theorem (Watson 1965]), from which we derive fi- 
nally 

F{a) ^ Io{ud/j'v)K (Lua/j'v) d<a (27) 
F(a) = Io(u)a/^'v)K (u!d/j'v) d> a (28) 
We deduce from (fTST) the (radial) electric field 

E(r,(j) 



27reon 2 u 2 7 



que 



iujz/v 



1 K {ujah'v)h(ud/ 1 'v) d<a(29) 
-floiwa/i'v) K^wd/Yv) d > a (30) 



E(r,uj) - , 

which has a discontinuity on the circle (as expected), of 
no physical consequence since it disappears on integration 
over the radius a. Comparing to the electric field spectrum 
(HHJ) of a point charge, we see that the finite size of the 
source removes the singularity of the field at d — > 0. It also 
increases the field strength outside the source by the factor 
Jo(wa/7'c), which increases as frequency increases, so that 
the lateral size of the charge makes the spectrum at large 
distances fall less quickly with frequency than that of a 
point charge. 
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Fig. 6. Spectral field strength | E(u) | at cj/2tt = 20 and 50 
MHz produced by 10 8 electrons with a NKG distribution 
with s = 1 (respectively 1.4) and tm — 100 m (respectively 
50 m), in solid lines (respectively dashed), for n = 1.0002 
and 7 = 40, as a function of perpendicular distance d to 
the charge's path, compared to the field of a point charge 
(thin lines). 



4.2. Pancake charge distribution 

Consider now a charge q of lateral distribution a (a) and 
finite longitudinal thickness. The potential at perpendicular 
distance d from the centre is obtained by integrating (|24p 

as 



/•OO 

$(r,w) = / da2iraa(a)G 
Jo 



(a) | $„(r,w) 



(31) 



for a distribution normalised so that 27T J °° daaa(a) = 1. 
Here G (a) involves an integral of the phase factor ^ z l v 
in Eq. ([24|) . ponderated by the longitudinal source distri- 
bution. The longitudinal extension therefore decreases the 
spectrum at frequencies uj > 1/St, where St is the width 
of the particle arrival time distribution, because of loss of 
coherence, so that we have G ~ 1 for u>5t < 1, whereas 
G <C 1 for ojSt ^> 1 when the fields no longer add in phase. 

We model the lateral distribution a{a) with the 
so-called NKG lateral density parametrization a (a) oc 

{a/r M ) s ~ 2 (1 + a/r M ) s " 4 - 5 (see Nishimura [HH7|) . Here the L 
two parameters s and tm are not necessarily the conven- 
tional shower age and Moliere radius respectively, since a 
wide set of couples [s, tm] can fit the observed lateral distri- 
bution (see Antoni & al. 120011 and refs. therein). However, 
the coupled variation of these parameters keeps roughly the 
same proportion of particles within a given lateral distance; 
for example, KASCADE measurements shown in Fig. 4 of 
the paper cited above find nearly 50% (respectively 80%) 
of the electrons within 20 m (respectively 50 m) lateral dis- 
tance. 

The distribution of arrival times is still less well known. 
Assuming 8t ~ 1.6 ns at the centre of the pancake, and a 
value roughly three times greater at 30 m lateral distance 
(Linslcy 1986), the factor G should decrease below unity 
at frequencies above about 50 MHz. In this highly simpli- 
fied calculation, we will therefore neglect this effect and the 
curvature of the pancake, keeping in mind that the corre- 
sponding loss of coherence should decrease the spectrum 
above about 50 MHz. 



Figure [5] shows the boosted Coulomb electric field spec- 
trum produced by 10 8 electrons with a NKG distribution 
with two sets of parameters: [s = 1,tm = 100 m] and 
[s = 1.4, tm = 50 m], at frequencies uj/2it = 20 and 50 
MHz, as a function of the perpendicular distance d to the 
charge's path, for n = 1.0002 and 7 = 40, compared to the 
field of a point charge. One sees that the degeneracy of the 
parameters of the lateral distribution only affects the field 
at very short radial distances, where it is smaller than that 
of a point charge by a factor depending on the parameters 
of the distribution. At distances greater than about 50 m, 
the field does not depend significantly on the details of the 
charge distribution, in the frame of this highly simplified 
model. 



4.3. Charge separation by the Lorentz force 

If there were no systematic separation of positive and 
negative charges by the Earth's magnetic field, the net 
charge producing the radio pulse would be equal to q ~ 
—eN(X)r]/2, where N(X) is given in (f2"3"]l and 77 is the rela- 
tive excess of electrons over positrons in the shower, which 
amounts to about 20 % at low energies. 

However, the Lorentz force due to the Earth's magnetic 
field B separates the positive and negative charges in the 
direction v x B. The boosted Coulomb electric field is thus 
the geometric sum of those produced by the positive and 
negative charges, each being directed along the perpendic- 
ular distance to their respective charge's barycentre, with 
opposite signs. The net electric field is then no longer di- 
rected along the perpendicular distance to the centre of the 
shower. For example, for a vertical shower impacting North 
of an antenna, this produces a Coulomb field having a East- 
West component, of strength roughly proportional to the 
separation between positive and negative charges. 

Let us estimate the average displacement L for electrons 
of Lorentz factor 7. The force F = —ev x B acting dur- 
ing the time t produces a displacement Ft 2 / (2jm e ) , in the 
plane perpendicular to v, directed normal to the projection 
of B in this plane. The time for travelling a radiation length 
Lq divided by the air density p is At = Lq/ (pc) , and the 
average of t 2 is 2 At 2 , so that the average displacement of 
electrons is 



eBL\ sin (3 
jrrieCp 2 



(32) 



where /3o is the angle between the magnetic field and the 
shower direction. With L ~ 367 kg/m 2 , so that L /(p c) ~ 
10~ 6 s at sea level, and B ~ 0.4 G, this yields 



L 



2100 sin/3 
l(p/Po) 2 



(33) 



With p = it/ 4 and 7 = 40, this yields L ~ 40 m at sea 
level and L ~ 60 m at an altitude of half the atmospheric 
scale height. Positrons are displaced in the opposite direc- 
tion, producing a separation of 2L. Allan (197TJ) finds values 
somewhat greater, but of the same order of magnitude. 

As an example, Fig. [7] shows the electric field spec- 
trum projected on the East- West direction, as a function 
of distance d to the centre along the NS axis, produced at 
u>/2tt = 10 and 50 MHz by 10 8 electrons and 10 8 positrons 
in vertical motion, displaced laterally by 50 m in opposite 
senses along the East- West direction, for n — 1.0002 and 
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Fig. 7. Spectral field strength projected on the East- West 
direction, produced at w/2ir = 10 and 50 MHz by 10 8 elec- 
trons and 10 8 positrons in vertical motion, displaced lat- 
erally by 50 m in opposite senses along the East- West di- 
rection, with NKG distributions with s = 1 and rj; = 100 
m (solid lines) compared to point charges (thin lines), for 
n = 1.0002 and 7 = 40, as a function of distance d to the 
centre along the NS axis. Exponential decreases are shown 
for comparison (dotted). 

7 — 40 (7' — 67); each charge species is assumed to have a 
NKG distribution with s = 1 and Tm = 100 m (solid lines). 
Exponential decreases are shown for comparison (dotted). 
Comparing to the case of point charges (thin lines), one 
sees that in this simple model, the field is not much af- 
fected by the NKG lateral distribution, except at very short 
distances. 

As an example of the change in polarisation introduced 
by the magnetic field, Fig. [8] shows contours of the East- 
West and North-South components of the boosted Coulomb 
electric field spectrum produced at 20 MHz by 10 s electrons 
and 10 s positrons in vertical motion, displaced laterally by 
50 m in opposite senses along the East- West direction; solid 
and dashed lines distinguish the sign of the field compo- 
nents; 7 and n have the same values as in FigJT) 

Note that the excess of electrons over positrons, not 
taken into account in this plot, produces an additional con- 
tribution changing both polarisation and amplitude. 

5. Electric field spectrum of a supraluminal charge 

Consider now a point charge moving faster than the waves, 
so that n/3 > 1. This will produce a Cerenkov field instead 
of a boosted Coulomb field. 



5.1. Retarded potentials 

As we noted in Section 2, the retarded potentials are still 
given by Eqs.fT])-©, but the geometry is different from 
Fig. [U since now PP' = nf3R > R. This means that P is 
on the other side of C on the z axis, at a distance so that 
cos > l/n/3. Therefore, the retarded source's distance R 
is now given by 



R 2 = d 2 + [n/3R - (z + vt)f 
which has two solutions if 
z + vt > d(n 2 j3 2 - l) 1/2 



(34) 



W field component 
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Fig. 8. Contour plots of the East- West and North-South 
components of the boosted Coulomb electric field spectrum, 
produced at w/27r = 20 MHz by 10 8 electrons and 10 8 
positrons in vertical motion, displaced laterally by 50 m in 
opposite senses along the Est- West direction, for n = 1.0002 
and 7 = 40. Contours levels correspond to strengths of 0, 

0. 1, 0.3, 1 (thicker lines), 3, 10 uV/m/MHz, with solid and 
dashed lines corresponding to components of different signs. 

1. e. t + z/v > r with r given by (|5|). When this inequal- 
ity holds, the antenna is inside the Cerenkov cone of half- 
angle arcsin (l/n/3) trailing the charge. When the oppo- 
site inequality holds, there is no solution for the retarded 
source, whereas when the two terms are equal (antenna on 
the Cerenkov cone), the two solutions merge together. 

This means that as time increases, the antenna first sees 
nothing (when it is still outside the Cerenkov cone), then 
it sees the field of one retarded source (when it is on the 
Cerenkov cone), which splits in two as the Cerenkov cone 
trailing the charge moves so that the antenna is inside it. 

Substituting into |TJ and summing on the retarded 
sources yields the potential 



*(r,t) 



2TTe n 2 d(n 2 [3 2 



a/2 



t + z/v 



for 

t + z/v > 



-1/2 
(36) 



(37) 



and 4>(r, t) = otherwise, which is similar to a result by 
Jackson (1999). The potential is singular at t + z/v = r, 
when the antenna is on the Cerenkov cone of the charge; 
at this time there is only one retarded source's position, 
satisfying z ret — z = d/ tan0 with cos 8 = l/n/3, i.e 



z ret -z = d (n 2 /3 2 - 1) 1/2 (38) 
so that the retarded altitude is now approximately 

2„2 i\-V2, 



h ret ~ d(n 2 (3 2 -iy 



cos a 



(39) 



for a shower of vertical inclination a with 7»l,n-l<l. 
We have plotted on Fig. [9] 



(35) 7c 



(n 2 /3 2 



1 



-1/2 



(40) 
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as a function of 



7 = (l — [3 Z ) " for different values of the refractive index 
under supraluminal conditions. 



as a function of 7. For n — 1 = 2 x 10~ 4 , particles with 
7 > 50 produce a Cerenkov field; however, one sees from 
Fig. |4] that, for the retarded altitude h ret ~ 7 c dcosa to 
ensure a significant retarded charge when d ~ 100 — 200 m, 
7 C should not be smaller than about 70; from Fig. [5] this 
requires 7 > 71 with the above value of n. 

This means that a little less than the higher half of the 
energy distribution can contribute to the Cerenkov field: 
the particles for which n(3 > 1 and whose retarded altitude 
is in the region of significant shower development. However, 
since the excess of electrons over positrons and the charge 
separation by the magnetic field ((33)) both decrease with 
energy, only a small proportion of these charges should con- 
tribute to the Cerenkov field. 

From (|36p . we deduce immediately the x component of 
the electric field E x = —d<&/dd 



E, r 



2ireon 2 d 2 (n 2 f3 2 



1 



■ V2 



t + z/v 



-3/2 



(41) 



when (|37p holds and when the opposite inequality holds. 



5.2. Cerenkov electric field spectrum 

Calculating the Fourier transform of the potential, which is 
given by (|36[) when the inequality (|3T|) holds and vanishes 
otherwise, yields the potentials in Fourier space 



$(r,w) 



iqe 



z/v 



4eon 2 i> 
A(r,w) = n 2 (3$(r,Lj)/c 



Hl{ud/ lc v) 



(42) 
(43) 



Here Hq = Jq + iYq is a Hankel function of the third kind 
and order l(Abramowitz & Stegun|1972), z is the antenna's 
co-ordinate along the charge's path (whose origin is the 
charge's co-ordinate at t = 0), and d is the antenna's per- 
pendicular distance to the charge's path. 

The electric field is given by substituting the potentials 
(|42 |) -(|43 |) into (fT5|) . so that its Fourier transform has the 
components 



iquje 



z/v 



4eon 2 v 2 "f c 
-que 1 "*/' 
4eon 2 v 2 "fl 



■Hl(wd/j c v) 
■Hl(uod/^ c v) 



(44) 
(45) 



where H{ = Ji + iY\. Therefore E z /E x < 1 for j c > 1, 
n — 1 <C 1, so that the electric field is nearly radial (from 
the charge's present position to the antenna), of amplitude 



iquje 



z/v 



(46) 



Aeon 2 v 2 "f c 

with 7 C given by 

Using the expansions of the Bessel function Hi 
(Abramowitz & Stegun ll972p . Eq. ([46|) yields at respectively 
low and high frequencies (or distances) 



cud , . qe kjz / v 

— « 1 E(r,u — 

7cf Zireon z vd 



(47) 



— » 1 E(v,co) 



qe 



i(uiz/v-K/i) /2ttujv\ 1/2 



2„2 



V led J 



for 7 C > 1, n - K 1. 

In the low- frequency (or small distance) limit, the elec- 
tric field has the same expression as for v < c/n, as expected 
since it is independent of 7 in this limit. This was already 
noted by Allan (|197ip . On the other hand, at large frequen- 
cies (or distances), the exponential decrease is replaced by 
a phase oscillation, again as expected since 7' is replaced 
by ij c when v > c/n. This phase oscillation suggests that 
the finite lateral size of the source should decrease the field 
with respect to that of a point charge, instead of increasing 
it as occurred for the boosted Coulomb field, because of loss 
of coherence. This effect is estimated in the next Section. 
Furthermore, the large frequency (or distance) value is es- 
sentially produced by the fast variation (see Eq.([H])) near 
t + z/v = t when the antenna is on the Cerenkov cone; this 
contribution is expected to be sensitive to the variation of 
the refractive index. 

6. Cerenkov spectrum of a charge distribution 

Consider a charge q distributed round a circle of radius a 
in the plane perpendicular to z which passes at z = at 
t = 0. The potential produced at perpendicular distance 
d from the centre of the circle (and distance z along the 
charge's path) is, from (|4"2"1) with the geometry defined in 
Fig(5] 

■ iuz/v 



4eon 2 i> 



1 



■F c {a) (48) 

(49) 

2ad cos (50) 

The integral F c (a) can be calculated analytically by ex- 
panding Hq in series of Bessel functions using Neumann's 
addition theorem (Watson 1966]), from which we derive fi- 
nally 



F c (a) = — I dcj>Hl(urh c v) 
2tt Jo 

,1/2 



F c (a) = Jo(ujd/^ c v) HQ(uja/^ c v) d<a 
F c (a) = Jo{Loa/^ c v) Hl{ujd/-) c v) d>a 



(51) 
(52) 



We deduce from (fT5|) the electric field along the perpen- 
dicular distance 



E x (r,u) 



iquje 



z/v 



4eo7i v 7, 

iquJ ^z/v 

Aeon 2 v 2 "f c 



■H l (uja/'-f c v)Ji(ujd/j c v) d<a (53) 
Jo(uja/j c v)Hl(u}d/j c v) d > a (54) 
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Fig. 10. Spectral field produced at frequency ui/2tt = 20 
MHz by 10 8 electrons with NKG distributions with s = 1 
and tm = 100 m (solid line) and s = 1.4 and Tm — 50 m 
(dashed line) for n — 1.0002 and 7 = 70 at a function of 
perpendicular distance d, compared to the field of a point 
charge (thin line). 



and along the charge's path 
qtde iuJZ / v 

E z (r,u>) = ^-^H (uia/^ c v)J Q (u>d/^ c v) d<a (55) 

quje iuJZ / v 

E z (t,u>) = 5-^-J (wa/7 c i;)i7 (wd/7 c w) d>a (56) 

The expressions (1551) - ([54")) are similar to those derived by 
Kahn & Lerche (1966) with a different formulation involv- 
ing the wave equation; note that the latter study did not 
include the boosted Coulomb field since the charge speed 
was assumed to be equal to that of the primary. 

As in the case v < c/n, the x component (radial in 
cylindrical co-ordinates around the charge's path) is much 
greater than the z one when j c 3> 1. 

Comparing to the electric field spectrum ([46]) of a point 
charge, we see that, as for the boosted Coulomb field, the 
finite size of the source removes the singularity at d — > 0, 
yielding E = at the centre (as it should by symmetry). 
Note also that, contrary to the Coulomb case, the finite 
lateral size of the source decreases the field strength, by 
the factor Jo(uia/j c v). 

This is illustrated in Fig. [10] , which shows the mod- 
ulus and phase of the Cerenkov electric field produced at 
frequency uj/2tt = 20 MHz by 10 s electrons with a NKG 
distribution with the two same sets of parameters as in 
Fig[51 for n = 1.0002 and 7 = 70 at perpendicular distance 
d = 100 m, compared to the field of a point charge (thin 
lines) . 

Another difference with the Coulomb field is that 
the phase varies rapidly with the charge energy, except 
at small frequencies and/or distances, since Hj(Q ~ 
(2/ttC) 1/2 e'tC-sV-i) for ( > 1; hence we see from (J54J) that 
integrating over a distribution of 7 should decrease the field 
strength. 

7. Comparison with the Cerenkov radiation 

We have seen that the electric field is essentially directed 
along the vector from the charge's present position to 
the antenna (x axis), i.e. perpendicular to the velocity. 



However, it is not this field that is at the origin of conven- 
tional Cerenkov radiation, but instead the (smaller) com- 
ponent parallel to the velocity (z axis). 

Indeed, the energy flow through a cylinder of radius d 
around the charge's path is given by 

dW -ldW e c f + °° , n ,_ „ 

— = — = — — / dz2ndE z B y (57) 

dz v dt v J_ oc 

since the wave magnetic field B — V x A with A given by 
(|4"5)) is along y. This integral, along the cylinder at a given 
time, may be transformed into an integral over time at a 
given point on the cylinder (using again dz = —vdt) 

dW f +oa 

— — = -2TTe cd / dtE z (t) B y (t) (58) 

J —oa 

Converting this time integral into a frequency integral, we 
obtain 

dW f + °° 

— = -2e cdReal / duiE z (uj) B* (uj) (59) 
dz j 

where B y (ui) = —dA z (uj) jdd. For a point charge q, $ (ui) 
is given by (l42|) and E z (uj) by (l45|) . Using the expansions 
for large arguments of the Hankel functions Hq and Hi , we 
obtain the power radiated by the particle 

(O.-.^.^p^,--^) ,60, 
dz AttcqC' 3 Jo \ n z v z J 

which is the well-known Frank- Tamm result. If instead of 
a point charge, the charge is distributed round a circle of 
radius a, we have from Eqs. ([35]). (|35]1. and J5S]) 

2 y r+00 / c 2 \ 

(P) = ^l ^(l-^)[Jo(ua M ? (61) 

with 7 C given by (|40|) - a result that can also be found 
directly using for example the general expression derived 
by Meyer- Vernet (|1988j) with a different formulation. 

The total electric field is thus much greater than would 
be suggested from an order of magnitude estimate using 
(frJTj) . Indeed (|6"Tj) involves the large distance values of E z 
and By = —vn 2 E x , whereas the antenna measures essen- 
tially a local value of E x , which is much greater than E z . 

8. Discussion 

8.1. Physical significance of the results 

We have calculated analytically the boosted Coulomb and 
Cerenkov contributions to the electric field spectrum pro- 
duced by an extensive cosmic ray shower, with a highly 
simplified model. These field contributions are due to the 
relativistic speed of the radiating charges and do not 
depend on the acceleration. The result may be under- 
stood more intuitively by using Feynman's formula E = 
—q/ (47reon 2 c 2 ) <i 2 i?/eft 2 xsign(l — n(3 ■ n), where 8 is the 
retarded angular position of the particle (Fig. []}. This for- 
mula, which can be derived from the Lienard-Wicchert po- 
tentials used in the present paper (Feynman [1964), was 
used by Allan (|1971[) in his seminal review. With Allan's 
notations (origin at t = when the charge passes at clos- 
est approach C), we have from Fig. [TJ t = dj (cs'm9/n) — 
dj (iitan^). The solutions for 9 (one solution for v < c/n, 
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two solutions in the Cerenkov cone for v > c/n), and 
the second derivatives are given respectively in Eqs.(A16)- 
(A17) and Fig. A5 of the paper by Allan (|1971|) . This yields 

an electric field varying as (l + t 2 /r 2 ) 3 ^ 2 for v < c/n and 

as -H(t - t) (< 2 /r 2 - l)~ 3/2 for v > c/n {H being the 
Heaviside function), in agreement with Eqs. (O and (|4"1"]) of 
the present paper. Feynman's formula is more intuitive than 
the Lienard-Wiechert formulation used here, but it leads to 
more complicated calculations for the present problem. 

Feynman's formula enlightens the fact that the boosted 
Coulomb and Cerenkov fields both stem from the strong 
acceleration of the retarded angular trajectory at relativis- 
tic speeds, even when the charge's speed is a constant. This 
takes place during a short interval r of time near the instant 
when the charge's path (not necessarily the charge itself as 
we already noted, see Feynman 1964) passes at closest ap- 
proach to the antenna (if v < c/n) or has the antenna on its 
Cerenkov cone (if v > c/n). As already noted, the boosted 
Coulomb and the Cerenkov field have the same spectrum 
at low frequencies (2wf < 1/r), which is the time integral 
of the field, and is independent on both the Lorentz factor 
7 and the refractive index n. 

8.2. Boosted Coulomb versus Cerenkov 

An important fact about cosmic ray showers is that the me- 
dian speed of secondary charged particles is roughly equal 
to the phase speed of radio waves in air at about half the 
atmosphere scale height. As a consequence, roughly the 
lower half of the particle energy distribution produces a 
boosted Coulomb field, roughly the higher half produces 
a Cerenkov field, whereas the particles moving nearly at 
the wave phase speed should not contribute except at very 
small distances, because in that case the retarded altitude 
of emission (respectively ((HJ and ([39]) for the Coulomb and 
Cerenkov cases) jumps above the atmosphere. 

The boosted Coulomb and Cerenkov field spectra of a 
point charge are equal at small distances and/or frequen- 
cies, but the Coulomb field decreases with frequency and 
distance more rapidly than the Cerenkov one. However, 
the Cerenkov field at large frequencies or distances is pro- 
duced when the antenna is very close to the Cerenkov 
cone, so that it should be very sensitive to the variation 
of the refractive index with altitude, much more so than 
the boosted Coulomb field which does not involve any sin- 
gularity. Furthermore, the Cerenkov field is produced by 
the high energy tail of the shower particles, for which the 
net charge and the magnetic separation are smaller. The 
Cerenkov field is thus expected to be smaller than the 
boosted Coulomb field - a question that should be examined 
more carefully. 

Finally, we have not considered the transition radiation 
(see Ginzburg & Tsytovich I1979|) that may be produced 
by the variation in refractive index, especially when the 
charges hit the ground. 

8.3. Comparison with synchrotron emission and with 
observations 

The present estimate does not take into account the electric 
field produced by the acceleration due to the Lorentz force 
(see Huege & al. 120071 and refs. therein), which yields op- 



posite (time varying) deviations of the velocities of charges 
of opposite signs. The field due to the acceleration by the 
Lorentz force has been calculated as a special case of syn- 
chrotron radiation (see Huege & Falcke I2005|) with each 
particle completing only a small fraction of gyration since 
the free path is much smaller than the gyroradius. 

However, it should be noted that the published simula- 
tions of that emission neglect the effect of refraction. Taking 
into account the air refractive index may change the spec- 
tral density and the radiation pattern, even for particles 
satisfying n(3 < 1, by introducing an "equivalent Lorentz 
factor" according to (fT0|) , whereas for nj3 > 1, the charges 
can catch up with the waves they emit. A naive applica- 
tion of the classical expression of the synchrotron formulae 
(Jackson 1999) would suggest that for a given charge, re- 
fraction should not change significantly the field at small 
frequencies or distances since in that case it does not de- 
pend on 7. However, refraction should change the retarded 
altitude of emission, and therefore the number of radiating 
charges, since it depends on the development of the shower. 
In particular, charges with n/3 ~ 1 have a retarded altitude 
above the atmosphere (except for very short perpendicular 
distances), and thus should not contribute to radiation. 

Nevertheless, it may be interesting to compare the 
boosted Coulomb field with the published simulations of 
synchrotron radiation. We find a Coulomb field strength 
(see for example Fig. [7]) that may not be negligible com- 
pared to the values calculated for synchrotron radiation (see 
Huege & al. 120071 and refs. therein). Furthermore it appears 
to have a spectral decrease with frequency and distance 
that is grossly similar; for example the scale of decrease 
is of the order of magnitude of 100 m at 20 MHz. More 
precise comparisons require further studies which are out- 
side the scope of this paper. Note however that the boosted 
Coulomb and the synchrotron spectral shapes have differ- 
ent physical origins. The roughly exponential decrease with 
both frequency and lateral distance is a generic property of 
the boosted Coulomb field, which stems from the expan- 
sion of the Bessel function K\ for large arguments (|2U1) . 
whereas the field strength tends to a finite constant (fT9|) 
as u> — » 0. In contrast, the synchrotron electric field spec- 
trum of a point charge does increase with frequency, with 
the field vanishing for u> — > 0, so that the field decrease 
with frequency exhibited by the synchrotron simulations is 
produced by the longitudinal extension of the source and is 
very sensitive to the time of arrival distribution (see Huege 
& Falcke UnnS]). 

The boosted Coulomb field strength shown in Fig.[7Jis of 
the same order of magnitude as the published observations 
(see Allanfl97Tl Ardouin [20061 Horneffer [20061 Lecacheux 
& al. I2007|) . including the scale of variation with distance. 
For example Fig. [71 calculated with a plausible source size 
and a total number of particles (roughly proportional to the 
primary energy) corresponding to shower maximum devel- 
opment for a primary of 2 x 10 17 eV, yields a low-frequency 
field Fourier transform of about 6 /iV/m/MHz at 100 m per- 
pendicular distance, decreasing with a spatial scale of about 
100 m; with the same parameters, the relative excess r\ of 
electrons over positrons yields an additional field of 77 times 
the value in Fig. [S) Note that the spectral densities (de- 
fined for positive frequencies only) are twice these values, 
and that the above estimates are only order-of-magnitude 
ones. Comparing with Fig. 38 of Huege & Falcke 120051 for 
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example suggests that this contribution may not be negli- 
gible. 

Finally, since the boosted Coulomb field is produced by 
the charges satisfying n/3 < 1 and whose retarded altitude 
(|SJ) (depending on n, too) is in the region of significant 
shower development, it depends on the refractive index. 
The field should thus be sensitive to atmospheric pressure, 
temperature, and humidity. 

8.4. Defects of the model 

In order to obtain analytical results, we made many sim- 
plifications which make our work little more than a pre- 
liminary step for more detailed calculations. In particular, 
one should take into account the variation of the refrac- 
tive index with altitude, the longitudinal distribution of 
the source, the variation of the charge distribution with 
(retarded) altitude, the deviations of the particle velocities, 
and the particle velocity distribution. Such studies, which 
require a detailed numerical simulation of the shower, are 
outside the scope of this paper. 

In particular, as indicated in Section 4, the longitudinal 
extent of the source should reduce the amplitude above 
about 50 kHz, because of loss of coherence. 

In order to obtain analytical results, we have consid- 
ered a fixed value of 7, close to the median value for shower 
secondaries. A realistic calculation should involve an inte- 
gration over the particle velocity distribution, taking care 
of the singularity at n/3 = 1 . This singularity is expected to 
be alleviated by the fact that the corresponding retarded 
altitude jumps above the atmosphere, so that the charges 
producing it cannot contribute. Hence, at distances of hun- 
dreds of metres and frequencies of tens of MHz, the inte- 
gration of the field strength over 7 is expected to produce 
a result behaving roughly as calculated in Sections 3 and 
4, with a value of 7' roughly equal to the maximum value 
for which the retarded altitude yields a significant charge. 

9. Conclusions 

We have evaluated the electric field produced by the rela- 
tivistic velocity of the charges in extensive cosmic ray show- 
ers, using a highly simplified model of the charge distribu- 
tion in order to obtain analytical results. We show that the 
refractive index plays an important role, not only in pro- 
ducing a Cerenkov-like field at supraluminal speeds, but 
also in changing the emission at subluminal speeds since 
the bulk of secondary particles move roughly at the wave 
phase speed. In particular, except for small perpendicular 
distances, this puts above the atmosphere the retarded alti- 
tude of the particles moving close to the wave phase speed, 
so that they do not contribute to the field. Even though our 
calculations are based on oversimplified hypotheses, the an- 
alytical formulas obtained are useful to estimate the impor- 
tance of various parameters and physical processes for the 
part of the field that does not depend on the acceleration. 
Our main results are listed below. 

1) A little less than the lower half of the charge energy dis- 
tribution in the shower produce a boosted Coulomb electric 
field. With plausible parameters of the charge distribution, 
this yields an electric field spectrum comparable to the val- 
ues calculated for synchrotron radiation and to those ob- 
served. In particular, we find that a primary of 10 17 eV 



should produce a spectral field of a few /LtV/m/MHz at 20 
MHz and 100 m perpendicular distance, decreasing roughly 
as an exponential with the product of perpendicular dis- 
tance by frequency, with a scale of decrease of about 100 m 
at 20 MHz. This is comparable to current measurements, 
and therefore does not appear to be negligible. 

2) Charges of higher energy produce a Cerenkov-like field. 
However, this field may be more difficult to detect for 
several reasons: first, the relative excess of electrons over 
positrons is small at these energies, making the net charge 
small; second, the magnetic separation between electrons 
and positrons decreases with energy as I/7; third, the phase 
varies with lateral position, so that the finite lateral size 
of the source decreases the field amplitude because of loss 
of coherence, except at very small frequencies; this is not 
the case of the Coulomb field, for which loss of coher- 
ence is only produced by the (small) longitudinal exten- 
sion. And finally, it is only at large frequencies or distances 
(u>d/^f' c v 3> 1) that the Cerenkov-like field may be greater 
than the boosted Coulomb contribution, and in this range 
the variation of the refractive index with altitude should 
play an important role. 

3) The radio electric field depends on the air refractive in- 
dex, and thus on the pressure, temperature and concen- 
tration of various components, especially humidity. This 
should introduce variations of radio emission with the alti- 
tude of emission and reception, and may produce significant 
day-night, seasonal and other effects. 

4) This study suggests several developments which are out- 
side the scope of this paper. In particular, taking into ac- 
count the variation of the refractive index with altitude, and 
making a detailed comparison with synchrotron emission. 
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